Control of the persistent currents in two interacting quantum rings through the 
Coulomb interaction and inter-ring tunneling 
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The persistent current in two vertically coupled quantum rings containing few electrons is studied. 
We find that the Coulomb interaction between the rings in the absence of tunneling affects the 
persistent current in each ring and the ground state configurations. Quantum tunneling between 
the rings alters significantly the ground state and the persistent current in the system. 



INTRODUCTION 

Mesoscopic metallic rings and nanoscopic quantum 
rings (QRs) are known for the Aharonov-Bohm effect 
and its persistent current [l| where the phase rigidity of 
the electron wave leads to quantum oscillations in the 
current. Recently, oscillatory persistent current was de- 
tected in self- assembled InAs/GaAs semiconductor QRs 
having a single electron^. Experimentally there have 
been a few approaches to fabricate coupled rings and 
little effort has been done on the system of vertically 
coupled QRsQ. Due to the ring geometry, these cou- 
pled ring complexes open a new route for measurement 
of quantum interference effects and for novel many-body 
states@, i, HQ. Two vertically coupled quantum rings 
(CQRs) form a new type of artificial molecule. The 
inter-ring distance and the tunneling strength together 
with the ring radius are new tunable parameters provid- 
ing more degrees of freedom to modulate and to con- 
trol the electronic structure and the persistent current 
of these ring shaped artificial molecules. In the present 
work, within the current spin-density functional theory 
(CSDFT)@,|, we study the persistent currents in the 
CQRs of few-electrons. The two coupled rings exhibit 
new molecular many-body states. The tunability of the 
inter-ring distance and the inter-ring tunneling strength 
leads to a very rich variety of many-electron ground states 
and persistent current spectra. The persistent current 
and magnetization in a single ring have been studied ex- 
tensively where the external magnetic field is the only 
parameter to alter the current at fixed ring radius. Here, 
we show that the persistent current in the CQRs depends 
on the inter-ring distance and the tunneling strength be- 
tween the two rings. Tunneling between the rings per- 
mits the exchange interaction between them inducing 
new quantum states that alter significantly the persistent 
current. This opens new ways to control the persistent 
current in quantum rings. 



THEORETICAL MODEL 

We firstly study the inter-ring Coulomb interaction ef- 
fects on the persistent current in two interacting QRs 



without tunneling. We model the QRs by a displaced 
parabolic potential V(r) — ^m*ujQ(r — r ) 2 in the xy- 
plane, where r = (x,y) = (r,9), loq is the confinement 
frequency and ro is the radius of the ring. The two 
stacked identical rings are separated by a distance d in 
the z direction. They interact with each other through 
the Coulomb potential. A homogeneous magnetic field 
B = Be z is applied perpendicularly to the xy-plane, 
which is described by the vector potential A = Bree/2 
taken in the symmetric gauge. The Kohn-Sham orbitals 
V-'jnmcr(r) = exp(— imO)(j) nrna {r)Z j (z) are used to express 
the density and ground state energy, where a =f or j is 
the z component of the electron spin and lo c = eB/m*c is 
the cyclotron frequency. The total density in the rings is 
P( r ) = J2a L™m l^nm CT (r)| 2 - In the z-direction this den- 
sity is approximated by a (5-function in each ring. The 
inter-ring Coulomb potential is given by 



V -inter 
H 



(r) = / dr' 



e 2 p(r')/2 



d|' 



(1) 



with the inter- ring distance d = |d|. The exchange- 
correlation potential is considered within the local- 
density approximation. 

In the CSDFT all the quantities are functionals de- 
pending on the spin-up (p') and spin-down (pJ-) densities, 
and the vorticity V(r) = x j p (r)/p(r). Therefore 

the exchange-correlation scalar and vector potential can 
also be written as a functional of these quantities [cl [Tol] . 
The paramagnetic current density j p (r) can be explicitly 
written as 
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The persistent current density, which is a measurable 
quantity, is given by j(r) = j p (r) + (e/m*c)p(r)A(r), 
where the second part corresponds to the diamagnetic 
current density. The persistent current is obtained by 
I = J j(r)dr = J p + /rf. We confirm that the magnetiza- 
tion in both the single and coupled QRs is proportional 
to the persistent current. 

When inter-ring quantum tunneling is permitted, the 
electronic states in different rings are coupled and the 
CQRs form an artificial molecule. There are bonding 
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and antibonding states with an energy splitting Asas 
between the two levels. Tunneling occurs through a po- 
tential barrier between the rings in the z-direction V(z) 
which is modeled by two coupled symmetric GaAs quan- 
tum wells of width 120 A and a barrier of 250 meV height 
between them. Within this model, the energy splitting is 
given by Asas = 22.86 exp[— 22.42d(ro)] meV as a func- 
tion of the inter-ring separation [4j. For the CQRs with 
inter-ring tunneling, each GS can be labeled by three 
quantum numbers (S z , M z , I z ): the total spin S z , total 
angular momentum M z and the isospin quantum number 
I,. 



NUMERICAL RESULTS 

In the calculations, we consider GaAs quantum rings 
with a ring radius ro = 2 ao and a confinement energy 
Tiuiq = 5 meV, where ao = y^h/m*ujo. The inter-ring 
Coulomb interaction effects are studied by varying the 
distance d between the two rings. The ground state (GS) 
of a single ring can be characterized by two quantum 
numbers: the total spin S z and the total angular mo- 
mentum M z . The persistent current is obtained for two 
interacting rings without tunneling between them, con- 
taining each 3 electrons, at different inter-ring distances. 
As shown in Fig. 1(a), the current is an oscillation func- 
tion of the magnetic field with semi-linear segments. The 
different segments appear because M z increases with in- 
creasing magnetic field leading to a GS configuration (or 
phase) transition as indicated in the figure. At small 
magnetic fields, the GS phase of a 3 electron QR is 
(S z , M 2 )=(l/2, 1). With increasing magnetic field, the 
total angular momentum increases as M 2 =l,2,4, and 5, 
but the total spin remains 5^ = 1/2. We observe that the 
position of the jumps in the current is not a monotonous 
function of d. Fig. 1(b) presents a phase diagram of the 
ground states of the QRs in the B vs. d plane. The hori- 
zontal dotted lines indicate the GS transitions of a single 
ring. The GS depends on the inter-ring distance d and 
this dependence is stronger when d is close to the diam- 
eter of the rings (2ro). It results mainly from the inter- 
ring Coulomb energy which is different for the different 
GS configurations (different M z ) as shown in the inset of 
Fig. 1(b). The energy difference A_E lntcr of the inter-ring 
coulomb energy -E s nt M between the different configura- 
tions (S Z ,M Z ) and (S' Z ,M' Z ) is plotted as a function of 
d: AE intcI = E^Z - E&t 1 ,, . The three curves are 
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zero for d ~ 2 ro. In comparison with the single ring, the 
phase transition occurs at lower (larger) magnetic fields 
when d < 2tq (d > 2ro). As a consequence, we found the 
remarkable behavior that the inter-ring Coulomb interac- 
tion affects significantly the persistent current. Fig. 1(c) 
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FIG. 1: (Color online) (a) The persistent current as a func- 
tion of magnetic field at different d in two interacting QRs 
containing each 3 electrons without inter-ring tunneling, (b) 
The phase diagram in the B-d plane. The dashed horizontal 
lines indicate the transitions in a single ring. Each phase is in- 
dicated by two quantum numbers (S Z ,M Z ). The inset shows 
the energy difference A_E mtcr of the inter-ring Coulomb en- 
ergy between the different configurations at a fixed magnetic 
field, (c) Variation of the persistent current in the rings as a 
function of d at B=0.7, 0.85, and 1.2 T. The horizontal lines 
are the current for the single ring case. 
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FIG. 2: (a) The total persistent current, (b) the paramagnetic 
and (c) the diamagnetic current as a function of magnetic field 
in the CQRs of d — O.lro and with N =3 electrons, (d) The 
same as (a) but now for N =6. The different phases are 
indicated by (S Z ,M Z ,I Z ). 

shows that, at B =0.7 and 1.2 T where the GS config- 
urations are the same for all d, the persistent current 
decreases with increasing d till a minimum and then in- 
creases slowly to its limit value of an isolated ring. How- 
ever, variation of the inter-ring distance at fixed magnetic 
field can also induce the GS phase transition and, conse- 
quently, a jump in the persistent current. This is the case 
of B = 0.85 T in Fig. 1(c) where the GS configuration 
changes from (1/2,2) to (1/2,1) at d =1.2r leading to an 
abrupt jump of the persistent current from -4 to 22 nA. 

The CQRs form an artificial molecule when inter-ring 
quantum tunneling is permitted. Fig. 2(a) shows the 
persistent current as a function of magnetic field in such 
an artificial molecule with a total of 3 electrons in the 
strong tunneling regime for d = 0.1 ro. In this regime, the 
CQRs behave as a single one because only the bonding 
states are occupied. This is evident by the isospin quan- 
tum number I z which is always 3/2. The spin and the 
angular momentum (S Z ,M Z ) follows the same sequence 
as those in Fig. 1(a). Figs. 2(b) and 2(c) give the cor- 
responding paramagnetic and diamagnetic currents, re- 
spectively. The paramagnetic current jumps with the 
GS phase change while the diamagnetic current is a lin- 
ear function of the magnetic field. The total momen- 
tum increases with increasing magnetic field with jumps 
AM Z = (N±l)/2. But the total spin and iso-spin remain 
S z = l/2 and I z = N/2. 

Fig. 2(d) presents the persistent current in the CQRs of 
6 electrons at d = 0.1 r which is in the strong tunneling 
regime where I z — N/2 — 3. When the electron number 
is even (iV=6), with increasing magnetic field, the to- 
tal angular momentum always increase by AM 2 = N/2 
when the system transits from one phase to the next as 
indicated in the figure. The total spin oscillates between 
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FIG. 3: The phase diagrams of the ground state configura- 
tions of the CQRs of N=6 electrons. The three quantum 
numbers (S z , M Z ,I Z ) indicate the ground state configuration. 



and 1. The most stable states correspond to the closed- 
shell configurations (S z = 0). When the total spin is 
non-zero (S z = 1), the ground state is less stable because 
it corresponds to an open-shell system. 

In the cases of a single ring and the strongly coupled 
double rings, new GS phases appear with increasing mag- 
netic field due mainly to the increase of the total angu- 
lar momentum. However, with increasing inter-ring dis- 
tance, when the tunneling splitting energy between the 
bonding and antibonding levels reduces to compete with 
the inter-ring exchange-correlation energy, the molecular 
phases can be complex and so the persistent current. 

In order to understand better the persistent current in 
the CQRs in the intermediate and weak tunneling regime, 
we plotted in Fig. 3 the phase diagram of the CQRs of 
6 electrons (on the average 3 electrons per ring) in the 
B — d plane. The phase diagram of the corresponding two 
ring case without tunneling was presented in Fig. 1(b). 
It is seen that for small d, the coupled rings behave like a 
single one because of the strong tunneling (large Asas)- 
All the electrons are in the bonding states in this case. 
With increasing d, the coupled rings are in the so-called 
molecular phase. For 0.15ro ;$ d < 0.4ro, competition 
between the tunneling energy Asas and the electron ex- 
change energy leads to a quite complex phase diagram. 
The electrons occupy the antibonding state but most of 
the electrons are in the bonding state (I z = 1 or 2). 
At large d, Asas approaches zero, the same number of 
electrons are in the bonding and antibonding states, so 
I z = 0. In this case, the inter-ring direct Coulomb inter- 
action and the exchange energy dominate the GS phases. 
In comparison with Fig. 1(b), the inter-ring tunneling al- 
ters the phase diagram dramatically. For instance, in the 
phase (3,6,0), 3 electrons are in the bonding and 3 in 
the antibonding state. All these 6 electrons have parallel 
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FIG. 4: (Color online) (a) The persistent current, (b) the 
paramagnetic and (c) the diamagnetic current in the CQRs 
with d = 0.23ro and N =3 electrons. The dotted (red) and 
dashed (blue) curves are currents carried by the bonding and 
antibonding states, respectively. The persistent current in the 
CQRs of N=6 electrons at (d) d = 0.23r and (e) d = 0.40r . 
The thin dash-dotted curves are the current in the corre- 
sponding coupled rings without tunneling. Each GS phase 
is indicated by (S z , M z , I z ) 



spin. In a real sample, the tunneling at large inter-ring 
separation can be controlled by a weak external electric 
field in the ^-direction. 

In the intermediate tunneling regime (e.g., d = 0.23 
i - o), the CQRs are in the molecular phase. Fig. 4(a) 
shows the persistent currents as a function of the mag- 
netic field in the CQRs of 3 electrons. In this regime 
both the bonding and antibonding states are occupied as 
indicated by the value of I z . Therefore, more combina- 
tions of angular momenta become possible as a function 
of the magnetic field leading to more oscillations in the 
persistent current. At large magnetic fields, the rapid 
variation of the angular momentum with increasing mag- 
netic field leads to a suppression of the amplitude of the 
persistent current. The dotted (red) and dashed (blue) 
curves indicate the currents carried by the electrons in 
the bonding and antibonding states, respectively. The 
black curve represents the total current, which is a sum of 
the currents due to the bonding and antibonding states. 



Figs. 4(b) and 4(c) show the corresponding paramagnetic 
and diamagnetic current, respectively. In Fig. 4(b), we 
see that the paramagnetic current due to the antibonding 
state is zero until B=2.35T where the (1/2,6,1/2) state 
appears. An increase of the current due to antibonding 
states is always partially compensated by a decrease of 
the current due to the bonding states. Such a compen- 
sation is complete for the diamagnetic current as shown 
in Fig. 4(c). 

Figs. 4(d) and 4(e) show the persistent current as a 
function of the magnetic field in the CQRs of 6 electrons 
at d = 0.23ro and 0.40ro, respectively. The GS phase di- 
agram of these CQRs is given in Fig. 3. The dotted-dash 
(green) curves in Figs. 4(d) and 4(e) give the current in 
the corresponding rings without tunneling (3 electrons 
in each ring). Clearly the tunneling leads to different 
ground state phases and consequently different persis- 
tent current. We notice that in the CQRs of d = 0.23r , 
the persistent current is zero at B = (see Fig. 4(d)) 
when there is tunneling between the two rings. But this 
current is about -20 nA when the tunneling is absent. We 
know that a vertical gate potential can control tunneling 
in such a structure. It means that an external gate can 
switch on/off the persistent current (and the magnetiza- 
tion) of the these coupled rings even in the absence of a 
magnetic field. 



CONCLUSIONS 

In summary, the only external parameter which con- 
trols the persistent current in a single QR is the external 
magnetic field. In this calculation we show that in cou- 
pled quantum rings, the persistent currents in the CQRs 
can be controlled by the inter-ring coupling. This inter- 
ring coupling consists of the inter-ring Coulomb inter- 
action and the inter-ring quantum tunneling. The for- 
mer depends not only on the inter-ring distance, but also 
on the electronic configuration of the QRs which can be 
modulated by an external magnetic field. The later de- 
pends on the inter-ring tunneling energy, the inter-ring 
electron exchange effects as well as the external magnetic 
field and can also be controlled by a vertical gate. 

This work was supported by FAPESP and CNPq 
(Brazil) and by the Flemish Science Foundation (FWO- 
VI) and the Belgium Science Policy (IAP), Belgian State 
Policy. 



* Electronic address: hai@ifsc.usp.br 
[1] Aharonov-Bohm and other cyclic phenomena, J. Hamil- 
ton (Springer- Verlag, Berlin, 1997). 
[2] N. A. J. M. Kleemans et ai, Phys. Rev. Lett. 99, 146808 
(2007). 



5 



[3] D. Granados, J. M. Garcia, T. Ben, and S. I. Molina, 
Appl. Phys. Lett. 86, 071918 (2005); E. Marega, Z. Y. 
Abuwaar, and G. Salamo, (to be published, 2007). 

[4] L. K. Castelano, G.-Q. Hai, B. Partoens, and F. M. 
Peeters, Phys. Rev. B 74, 045313 (2006), and references 
therein; G. Piacente and G.-Q. Hai, J. Appl. Phys. 101, 
124308 (2007). 

[5] F. Malet, M. Barranco, E. Lipparini, R. Mayol, M. Pi, 
J. I. Climente, and J. Planelles, Phys. Rev. B 73, 245324 
(2006); M. Royo, F. Male t, M. Barra nco, M. Pi, and J. 
Planelles, Cond-Mat larXiv:0809.3"T03l (2008). 



[6] B. Szafran, Phys. Rev. B 77, 235314 (2008); Phys. Rev. 

B 77, 205313 (2008). 
[7] Y. Z. He and C. G. Bao, J. Phys.: Condensed Matt. 20, 

055214 (2008). 

[8] G. Vignale and M. Rasolt, Phys. Rev. Lett. 59, 2360 

(1987); Phys. Rev. B 37, 10685 (1988). 
[9] S. Viefers, P. Singha Deo, S. M. Reimann, M. Manninen, 

and M. Koskinen, Phys. Rev. B 62, 10668 (2000). 
[10] B. Partoens and F. M. Peeters, Europhys. Lett. 56, 86 

(2001); Phys. Rev. Lett. 84, 4433 (2000). 



